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ABSTRACT 

We critically re-examine the available data on the spectral types, masses and radii 
of the secondary stars in cataclysmic variables (CVs) and l ow-mass X-ray binaries 
(LMXBs), using the new catalogue of Ritter & Kolb (1998) as a starting point. We 
find there are 55 reliable spectral type determinations and only 14 reliable mass de¬ 
terminations of CV secondary stars (10 and 5, respectively, in the case of LMXBs). 
We derive new spectral type-period, mass-radius, mass-period and radius-period re¬ 
lations, and compare them with theoretical predictions. We find that CV secondary 
stars with orbital periods shorter than 7-8 hours are, as a group, indistinguishable 
from main sequence stars in detached binaries. We find it is not valid, however, to 
estimate the mass from the spectral type of the secondary star in CVs or LMXBs. We 
find that LMXB secondary stars show some evidence for evolution, with secondary 
stars which are slightly too large for their mass. We show how the masses and radii of 
the secondary stars in CVs can be used to test the validity of the disrupted magnetic 
braking model of CV evolution, but we find that the currently available data are not 
sufficiently accurate or numerous to allow such an analysis. As well as considering sec¬ 
ondary star masses, we also discuss the masses of the white dwarfs in CVs, and find 
mean values of M i = 0.69 ± 0.13Mq below the period gap, and M\ = 0.80 ± 0.22M@ 
above the period gap. 


Key words: binaries: close - novae, cataclysmic variables - X-rays: stars - late type 
stars. 


1 INTRODUCTION 


Catacl ysmic variables (CVs) arc gcmi-dctachcd binary stars 
in which a white dwarf primary accretes material from a 
Ro che-lobe filling se condary. For a thorough review of CVs 
see Warner (1995a). The spectral type and luminosity class 


of the secondary star can be estimated from basic theory, as 
follows. Kepler’s third law can be written as 


4-7r 2 a 3 , . „ , ( 1 + q\ 

_ F =M 1+ M2 = M 2 (—- J, 


(1) 


where q = M 2 /Mi. The volume-equivalent radius of the 
Roche lobe can be approximated by 


R 2 

a 


-At hC 


( 2 ) 


This relation is a slightly modified form of the Paczyriski 


(1971) equation, and we have found that it is accurate to 


less than 3 per cent over the range of mass ratios relevant 
for CVs (0.01 < q < 1.0). Combining equations (1) and (2) 
gives the mean density-period relation, 


PQ 


/ m 2 \ 

( R 2\ 

\Mq) 

\r q J 


= 75.5P“ 2 (hr), 


(3) 


which is accurate to ~ 6 per cent (see Eggleton 1981). 

Typical lower main-s equence me an densities range from 
~ 5Op0 for M8 dwarfs (Allen 1976), corresponding to the 
minimum orbital period of CVs around 80 min, to ~ 1 pg 
for GO dwarfs, corresponding to an orbital period of around 
9 hr. At longer periods than ~ 9 hr, the density of the sec¬ 
ondary star must be sub-solar, corresponding to F-type (and 
earlier) main sequence stars. F-type main sequence stars 
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have masses above solar (M ~ 1.3Mq for an F5 dwarf, _ properties. A similar conclusion has been reached by Ritter 


M ~ 1.6M© for an FO dwarf, Gray 1992), which is rele¬ 
vant when considering the formal requirem ent of q < 5/6 
for sta ble, c onservative mass transfer (e.g. Frank, King & 
Raine 1992). Since the mass of the white dwarf must be 


below the Chandrasekhar limit, the secondary is forced to 
have M2 < 1.2 Mq. As the white dwarf population is biased 
towards much lower masses (the mean white dwarf mass in 
CVs is 0.77±0.21Mq, see section 6), the number of systems 
with secondary stars above solar mass is very small indeed. 
Secondary stars with sub-solar mean densities in long pe¬ 
riod CVs must therefore be low mass evolved M or K-stars 
rather than intermediate mass main sequence F-stars. In 
summa ry, for the orbital period range in which most CVs 


lie (1.3 rr<P<9hr), the secondary stars should be M, I< or 
G main-sequence dwarfs, while longer period systems must 
harbou r secondaries which have evolved away from the main 


sequent e. 


Even though we have just shown that most CV secon¬ 
daries should have lower main-sequence densities, it is not 
clear that they should appear as lower main-sequence stars. 
This is because CV secondaries are subject to a number of 
extreme environmental factors to which held stars are not. 
Specifically, CV secondaries are: 

17?© from a hot, irradiating source (see 


1-11 AT© yr -1 , 


(i) situat ed 
Smith 199f| ), 

(ii) rapidly rotating (~ 100 kms -1 ) 

(iii) Roche-lobe shaped, 

(iv) losing mass at a rate of ~ ICE 8 — 10“ 

(v) survivors of a common-envelope phase during which 
they existed within the atmosphere of a giant star, 

(vi) exposed to nova outbursts every ~ 10 4 yr. 

It is the purpose of this paper to see if the above envi¬ 
ronmental factors alter the gross properties (masses, radii, 
spectral types) of the secondary stars in CVs. This is not 
the first time that the question of whether CV secondaries 
lie on t he main sequence has been investigated. Echevarria 


(1983) concluded that they were not. He came to this con¬ 


clusion by calculating an empirical mass-radius relation for 
held stars and then using a density-period relation similar 
to equation (3) to calculate a mass-period relation. This was 
then combined with a spectral type-mass relation to get a 
spectral type-period relation for held stars. By comparing 
this relation with the measured spectral types and periods 
of 17 CVs he found that the CVs did not Ht his relation 
and concluded that they were therefore not on the main 
sequence. 

Warner (1995a,b) disputed Echevarria’s conclusion by 
arguing that he had equated a poorly-htting power law to 
the mass-spectral type relationship and that as a result 
of this Echevarria’s derived relationship systematically pre¬ 
dicted spectral types that were too early for long period 
CVs and too late for the short period CVs. Instead, Warner 
plotted the CV secondary stars on a spectral type-period di¬ 
agram and calculated equivalent periods of held stars using 
a mass-period relation, derived us ing a mass-radi us rela¬ 
tion which hts b oth the CV dat a of Webbink (1990) and the 
held star data of Popper (1980)| . Warner concluded that, as a 
group, secondary stars have masses, radii and spectral types 
related in exactly the same way as main sequence stars, al¬ 
though some individual CVs do depart from the average 


(1983) and Patterson (1984) 


Low mass X-ray binaries (LMXBs) are also semi¬ 
detached systems, the difference being that the primary in 
these systems is a neutron star or a black hole. The sec¬ 
ondary stars in LMXBs should be similar to those in CVs, 
except that they are irradiated to a much greater degree by 
X-rays. Also, the formation of LMXBs containing neutron 
stars requires that the secondary stars have a mass between 
1.3 — 1.5 Mq at the onset of mass tra nsfer, and that the y 
should be significantly nuclear evolved ( King fe Kolb 1997 ). 
LMXB secondaries should therefore appear even less like 

main sequence stars than their CV counterparts. _ 

In this paper we repeat the analy ses of [Echevarria 


(1983) , Ritter (1983) , Patterson (1984) , Webbink (1990) 


and Warner (1995a,b) and go a few steps further, using 
the m ore extensive data set now available (Ritter & Kolb 


1998 and other sources) to derive the first reliable spectral 


type-period relation for CVs. We present new mass-radius, 
mass-period and radius-period relations for CV secondary 
stars with reliable system parameters, and compare the pa¬ 
rameters of the CV and LMXB secondaries with those of de¬ 


tach ed binaries as l isted in the reviews of Andersen (1991) 


and Popper (1980). We also briefly examine the observed 
mass distribution of the white dwarf population in CVs, and 
conclude with a discussion on how the measured masses and 
radii of CV secondary stars can be used to constrain the dis¬ 
rupted magnetic braking model of CV evolution. 


2 THE ORBITAL PERIOD-SPECTRAL TYPE 
RELATION 

The spectral type of the secondary star has been measured 
in 66 CVs with known orbital per iods, according to the cat¬ 


alogue of Ritter & Kolb (1998) A number of these mea¬ 


surements, however, are dubious, relying on main sequence 
assumptions, o r infrared colours which are contam inated by 
disc emission (Berriman, Szkody & Capps 1985), to deter¬ 
mine the spectral types. The only reliable method of de¬ 
termining spectral types is to detect absorption features in 
the spectra of the secondary stars. One can then compare 
the equivalent widths of the secondary star absorption fea¬ 
tures with those of isolated dwarf stars (e.g. Wade & Horne 
1988), or simply match the observed absorption features 


with those in a set of template spectra, either by eye (e.g. 
Friend 199 0a,b) or by using an optimal subtraction tech¬ 
nique (e.g. [Smith, Dhillon fe Marsh 1998 ). Absorption fea¬ 
tures have been detected (and hence reliable spectral types 
determined) in 55 of the 66 CVs catalogued by Ritter & Kolb 


(1998), as listed in Table 1. Most of the spectral types listed 


in Table 1 do not have error bars associated with them, so 
we have instead assigned a weighting to each determination. 
In most cases, a weighting of 1 implies a careful spectral type 
determination using ample template dwarfs and a weighting 
of 0.5 implies there is a larger uncertainty in the exact value 
of the spectral type, often as a result of insufficient template 
dwarfs having been observed. 

The spectral types of the secondary stars listed in Ta¬ 
ble 1 are plotted against orbital period in Fig. 1. Two linear 
least-squares hts were performed, one above and one below 
the kink at an orbital period of 4 hours. Note that we omit- 
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Orbital Period (hr) 

Figure 1. The spectral types of the secondary stars in CVs and LMXBs versus their orbital periods. The CVs have been plotted as 
filled circles, large circles denoting better spectral type determinations (weight = 1 in Table 1) than the small circles (weight = 0.5). 
LMXBs are denoted by triangles. Also plotted are 50 isolated stars for which the spectral types and masses have been measured (see 
text for details). The solid line is a two part linear fit to the CVs, given by equation (4). The dashed line is the relationship derived by 
Echevarria (1983). 


ted 7 CVs from the fit (indicated by the asterisks in Table 1) 
which are believed to harbour evolved secondary stars. For 
the sake of simplicity we have treated each spectral type as 
being divided into ten equal sub-types.Q The resulting fits 
are given by: 


26.5 

± 0.7 

- 0.7 P, 
± 0.2 

P < 4hr 

33.2 
± 3.1 

- 2.5 P, 
± 0.5 

P > 4hr 


the dashed line represents Echevarria’s equation. Also plot¬ 
ted as a comparison in Fig. 1 are the Sun and 49 other 
late-type main sequence stars in detached binaries. Twelve 


of these 49 have been taken from the list of 

Andersen (1991) 

while the review paper by 

Popper (1980) 

lists a further 37 


stars with less well refined mass and radius measurements, 
but which are suitable for comparison purposes; we have 
taken 11 stars from Popper’s Table 2 (Detached main se¬ 
quence binaries, B6 to M), 6 stars from Table 7 (Resolved 
spectroscopic binaries), and 20 from Table 8 (Visual bina¬ 
ries) . 


where Sp( 2) = 0 represents a spectral type GO, Sp( 2) = 10 
is K0 and Sp( 2) = 20 is M0. The rms scatter is 0.8 spectral 
sub-types for P < 4 hr, and 3.0 sub-types for P > 4 hr. 

The solid lines in the spectral type-period diagram rep¬ 
resent our two-part linear fit to the CVs (equation 4), while 


* The MK system officially does not contain the K8 or K9 sub - 
types, nor those of M4.5, M5.5 etc. (Jaschek & Jaschek 1987). 
However these are in (infrequent) use, so to make the spectral 
type—period relation simple to use, and since the system is in any 
case non-linear, we have adopted a system where each spectral 
tv nc is divided into 1 0 si ihtvnos. This sys tem was also followed 
by Echevarria (1983) and Warner (1995b) in their spectral type- 
period diagrams. 


The stars in detached binaries are plotted in Fig. 1, 
not according to their actual orbital periods, but rather ac¬ 
cording to the period of a CV containing a Roche lobe-filling 
secondary star of that mass and radius. This assumes, there¬ 
fore, that the spectral types of the stars in detached binaries 
would remain unchanged if they were to become lobe-filling 
secondaries. There are two ways of performing the conver¬ 
sion from mass and radius to period, depending on how ac¬ 
curately the radii of the detached stars are known. In the 
case of the stars listed by Popper, which have poorly deter¬ 
mined radii, we have converted from mass to period using 
the empirical mass-period relation we derive in section 4 
(equation 8). Typical errors in this conversion, as calculated 
from the errors in equation (8), are plotted for different or¬ 
bital periods at the top of Fig. 1. In the case of the stars 
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listed by Andersen, all of which are eclipsing binaries and 
hence have accurately determined radii, we have used equa¬ 
tion (3) directly to convert from mass and radius to period. 
These have a much smaller associated error (3 per cent), 
so the crosses in Fig. 1 (Andersen’s points) should be given 
greater weighting than the open circles (Popper’s points) 
when interpreting the diagram. 

There appears to be little difference between the distri¬ 
butions of the CVs and the main sequence stars in Fig. 1 
for periods below 4 hours. Above 4 hours, there is much 
larger scatter in the two distributions but they still seem to 
agree up to 7-8 hours. Beyond 7-8 hours the two distribu¬ 
tions show a definite divergence and some evolved secondary 
stars appear, while beyond 9 hours the secondary stars all 
show signs of evolution (as is to be expected, see section 1 ). 
Echevarria’s relationship is a poor fit to the detached sys¬ 
tems, as noted by Warner (1995a,b), and can be seen to er¬ 
roneously predict later spectral types at short periods, and 
earlier spectral types at long periods, than those which have 
been observed. Note that the Echevarria relationship was 
used as the basis for the claim by Friend et al. (1990a,b) 
that the secondary stars in CVs were “too cool for comfort” 
and “too cool for credibility,” which no longer appears to be 
true (but see Section 4). 

In Fig. 1, we have also plotted the 5 low-mass X-ray 
binaries (LMXBs) with known orbital periods (below 12 
hours) and spectral types. They are listed in Table 2. To 
within the scatter, the LMXBs appear to have the same dis¬ 
tribution as the CVs, indicating that LMXBs with orbital 
periods of 7-8 hours and longer harbour evolved secondary 
stars. 


3 THE MASS-RADIUS RELATION 


The masses of the component stars in CVs can be deter¬ 
mined in a number of different ways, using measurements 
of: 

(i) The radial velocity semi-amplitude of the white dwarf, 
Kw- This is inferred from the velocity variations of the wings 
of the emission lines, which arise in the inner accretion re¬ 
gions and are assumed to follow the motion of the white 
dwarf. These measurements are often unreliable because of 
contamination from other emission-line sources, such as the 
bright spot and the irradiated face of the secondary star, 
which do not follow the motion of the white dwarf and hence 
introduce phase shifts in the radial velocity curves. It is pos¬ 
sible to correct for this contamination using diagnostic di¬ 
agrams ( shatter, Szkody fe Thorstensen 1986 ), light centre 
diagrams ( Marsh 1988 ) and symmetry analyses of Doppler 
tomograms ( ptilfT996 ). However we do not accept that any 
of these methods do give truly reliable values of Kw because 
of the uncertainty in the extrapolation needed to correct for 
the large phase shifts in the radial velocity curves. 

(ii) The radial velocity semi-amplitude of the secondary 
star, Kr. This is measured using the motion of absorp¬ 
tion lines such as the Nal A8190A doub let (e.g. Friend et 
al. 1990a,b) or by using skew mapping (Smith, Dhillon & 


Marsh 1998 ). Another option is to follow the motion of the 
secondary sta r using the line emission from its irradiated in¬ 
ner face (e.g. Beuermann & Thomas 199C). All of these Kr 
measurements are subject to errors due to the non-uniform 


distribution of the emission/absorption line strength on the 
surface of the secondary star, although this can be cor rected 
for (e.g. Wade & Horne 1986, Rutten fe Dhillon 199^ ). 

(iii) The projected rotational velocity of the secondary 
star, V sin i. This is usually measured through comparison 
with (slowly rotating) field star template spectra which are 
given a series of artificial rotational broadenings. A constant 
times each template spectrum is then subtracted from the 
object spectrum to give a residual spectrum. The template 
which yields the smoothest residual provides the value of 


V sin i (e.g. Marsh, Robinson fe Wood 1994 


Smith, Dhillon 


fc Marsh 1998[ ). An alternative method of measuring V sin* 
is that of cross-correlating with a template spectrum and 
measuring the width of the c ross-correlation peak, a ro t (e.g. 
Horne, Wade & Szkody 1986). Models of rotating stars with 
various amounts of limb darkening are then used to find val¬ 
ues of V sin* corresponding to a ro t . Measurement of V sin* 
provides a powerful constraint on M 2 , as the relation be¬ 
tween them is only weakly dependent on q (see equation 19 
in Appendix B). 

(iv) The full width of the eclipse at half depth A^j/ 2 - 
At half depth it is assumed that one half of the accretion 
disc is eclipsed, corresponding to the point at which the 
white dwarf is eclipsed. This ceases to be true if the disc is 
asymmetrical, e.g. if there is a dominant bright spot, or if 
there is no disc, e.g. in magnetic CVs. If the system is non¬ 
magnetic and does not have a dominant bright spot then 
the measurement of A 0 J /2 provides a relationship betwee n 
q and the inclination i (e.g. Smith, Dhillon & Marsh 1998). 

(v) The duration and phase of the ingress/egress of the 
bright spot. This assumes that the bright spot lies on a bal¬ 
listic trajectory from the L\ point, i.e. there is no significant 
magnetic fie ld. This provides another relationship between 
q and i (e.g. Wood et al. 1989). 

(vi) The duration of the ingress/egress of the white dwarf, 
A (j> w d- This is a measure of the radius of the primary as 
a function of i and the separation a. The presence of an 
extended boundary layer can, however, cause the radius of 
the white dwarf to be ov er-estimated and h ence the mass to 
be under-estimated (e.g. Wood et al. 1989 ). 

(vii) The radius of the orbit of the white dwarf about the 
centre of mass, a w d, for rapidly rotating white dwarfs. This 
can be estimated from the spin pulse delay (e.g. AE Aqr, 
Eracleous et al. 1994 ). The only question is whether the 
white dwarf is the source of the pulses; in AE Aqr the spin 
pulse delay is exactly in phase with the white dwarf so there 
is little doubt. 

(viii) The ellipsoidal variations due to the changing as¬ 
pect of the distorted Roche lobe-filling secondary star. These 
are particularly prominent in the infrared, showing up as a 
distinctive double-humped modulation in the light curve. 
With appropriate modelling it is possible to constrain the 
inclination using ellipsoidal variations. This technique is es - 
pecially useful for non-eclipsing systems (e.g. Hilditch 1995). 
The main uncertainty with this technique is the contribution 
of the disc to the total flux, which reduces the amplitude of 
the variations, resulting in an underestimate of the inclina¬ 
tion. This can be corrected for i f the secondary star ha s been 
detected spectroscopically (e.g. Sliahbaz et al. 1996a). 

(ix) Linear polarisation light curves o f magnetic CVs. 
These can be modelled to constrain i (e.g. Wickramasinghe 
et al. 1991 ). 
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Table 1. Spectral types of secondary stars and orbital periods (in hours) of cataclysmic variables. Key: CN = classical nova; DN = 
dwarf nova; RN = recurrent nova; NL = nova-like; P = polar; IP = intermediate polar; * = systems with ” evolved” secondary stars, 
omitted from the linear fit in Fig. 1. 


Star 

Class 

Porb 

Sp Type 

Weight 

Ref. 

Star 

Class 

^orb 

Sp Type 

Weight 

Ref. 

BZ UMa 

DN 

1.63 

M5.5V 

1 

4,5 

DQ Her 

CN,IP 

4.65 

M3+V 

1 

42,43 

RX J0719+655 

P 

1.63 

M4-6V 

0.5 

1 

UX UMa 

NL 

4.72 

MOV 

0.5 

44,45 

EX Hya 

IP 

1.64 

M3V 

0.5 

6,7 

V895 Cen 

P 

4.77 

M2V 

1 

74 

V834 Cen 

P 

1.69 

M6.5V 

0.5 

8,9 

EX Dra 

DN 

5.04 

M1-2V 

1 

46,47 

HT Cas 

DN 

1.77 

M5.4V 

1 

10,11 

RX And 

DN 

5.04 

K5V 

0.5 

48,49 

Z Cha 

DN 

1.79 

M5.5V 

1 

12,13 

AR Cnc 

DN 

5.15 

M4-5.5V 

0.5 

2,3 

V2301 Oph 

P 

1.88 

M6V 

1 

14,15 

EY Cyg 

DN 

5.24 

K5-M0V 

1 

50,28 

MR Ser 

P 

1.89 

M5-6V 

0.5 

16,17 

CZ Ori 

DN 

5.25 

M2.5V 

1 

4 

ST LMi 

P 

1.89 

M5-6V 

0.5 

17,18 

AT Cnc 

DN 

5.73 

K7-M0V 

1 

51,28 

AR UMa 

P 

1.93 

M6V 

1 

19 

AH Eri 

DN 

5.74 

M3-5V 

1 

72,73 

DV UMa 

DN 

2.06 

M4-5V 

1 

20,3 

AH Her 

DN 

6.20 

K5V 

0.5 

52,53 

HU Aqr 

P 

2.08 

M4V 

1 

21 

SS Cyg 

DN 

6.60 

K5V 

1 

54 

QS Tel 

P 

2.33 

M4.5V 

1 

22 

V426 Oph 

DN 

6.85 

K3V 

1 

55 

AM Her 

P 

3.09 

M4+V 

1 

23 

Z Cam 

DN 

6.98 

K7V 

0.5 

56,57 

MV Lyr 

NL 

3.20 

M5V 

1 

24,25 

EM Cyg 

DN 

6.98 

K5V 

1 

58,59 

V1432 Aql 

P 

3.37 

M4V 

1 

26 

AC Cnc 

NL 

7.21 

K0V 

0.5 

60,61 

UU Aql 

DN 

3.37 

M2-4V 

1 

27,28 

TT Crt 

DN 

7.30 

K5-M0V 

0.5 

62 

QQ Vul 

P 

3.71 

M2-4V 

0.5 

29,17 

V363 Aur 2 

NL 

7.71 

G8V 

0.5 

63 

IP Peg 

DN 

3.80 

M4.5V 

0.5 

30,31 

V1309 Ori* 

P 

7.98 

M0-IV 

- 

64,65 

VY For 

P 

3.80 

M4.5V 

1 

32 

BT Mon 

CN 

8.01 

G8V 

0.5 

66 

KT Per 

DN 

3.90 

M3.3V 

1 

33 

CH UMa* 

DN 

8.23 

MOV 

- 

54 

CN Ori 

DN 

3.92 

M4+V 

1 

34,35 

RU Peg 

DN 

8.99 

K3V 

1 

67,54 

DO Dra 

DN,IP 

3.97 

M4V 

1 

36,80 

AE Aqr* 

IP 

9.88 

K4V 

- 

68,69 

WW Cet 

DN 

4.22 

M2.5V 

1 

37,38 

DX And* 

DN 

10.6 

K1V 

- 

70 

U Gem 

DN 

4.25 

M4+V 

1 

39,35 

U Sco* 

RN 

29.5 

F8V 

- 

75,76 

BD Pav 

DN 

4.31 

K7V 

1 

35 

GK Per* 

CN,IP,DN 

47.9 

K2-3IV-V 

- 

77,78 

TW Vir 1 

DN 

4.38 

M5-6V 

0.5 

40,28 

V1017 Sgr* 

CN,DN 

137.1 

G5IIIp 

- 

79,71 

SS Aur 

DN 

4.39 

M1V 

1 

41,35 








References: 1. 

5. 

10 . 

14. 


Ibvmassian et al. 1997. 2. Howell et al. 199o| , 3. Mukai et al. 1990 , 4. Ringwald, Thorstensen &; Hamwey 1994 


Sterken et al. 1983. 7. [D hill on et al. 1997. 8. Schw one et al. 1993 ~ 9. [Puclm arewicz et al. 1990 
4ornc. Wood & Stcining 199l[ 11. [Marsh 199p[ 12. [Robinson et al. 199fj. 13. Wade fc Horne 


Turcevic et al. 1 


42. Zhang et al. 1995. 43. Young fc Schneider 1981. 44. Baptista et al. 1995[ 45 


47. IBillington, 
51. 


Marsh fc Dhillon 1996 


25. Schneider. Young fc Scliectman 198l|. 26. [Watson et al. 1995[ 27. fitter fc Kolb 199£ , 2; 

29. Andronov fc Fuhrmann 1987|. 30. Wolf et al. 1993 31. ^'artin. .Tones fc Smith 1 )87, 32 
33. Thorstensen fc Ringwald 1997|. 34. Barrera fc Vogt 1989 , 35. [Friend et al. 1990a. 36. |H~asw ell et al 




*arwi^ 


Ritterj^Baernbantner^gg j, 15. Silber et al. 19941 16. gchwope' et a l~l_99l[ 17. [Mukai fc Charles 1986 

1 Te 1 


18. Cropper 1986|. 19. |R emi na^^^chachter fc Silbe^H194r20. [Ho well et, ai^K)88p21. |Cllenn~e^a^^fi94 . _ 

22. j^chwooe et al. 199f. 23. [Young. Schneider fc Sehectman 198 1| . 24. Skillman. Patterson fc Thorstensen 199. 


Smith et al. 199 1 


i 


euermann et a l. 1989 

nn 


37. linpwald of. ah 1 HKI 38. [HawkiiiS. Siiiit.li ,V Jones 1 99(7 39. Srna.k 199.^ 40. Sha.ft.or 1983 -II. Jhafter b. Harkncss 1986 

fiuttcn ct al. 199l| 46. [Fiedler 199l| . 


62. Jzkodv ct al. 199S. 63 


48. [Kaitchuck 1989[ 49. frlhilloii fc Marsh 199s|. 50. Jarna. Pvch fc Smith 1995 

_ .oetz 198f. 52. Horne. Wade fe Szkodv~D986[ 53. [Bruch 19871. 54. Friend ct al. 1990b[ 55. [kessman 198f| , 

56. rhorstensen V Ringwald 199fl . 57. |Szkodv fc Wade 1981. 58. [Stover. Robinson fc Nathcr 198ll _ 

59. 3euermann fc Pakull 1983. 60. Okazaki. Kitamura & Yamasaki 1982. 61. [Schleae]. Kaitchuck fc Honeycutt 1984 


Schleg el. H oneycutt fc K aitc huck 198f. 64. Buckley fc Sha tter 


66. jmith. Dhillon fc Marsh 199?, 67. 


70. Drew. Jones fc Woods 1993. 71. jckiguchi 199: 


Stover 1981 


3 68 . ^ 

2_72. pT 


r elsh. Horne fc Gomer 199: 


74. jtobie et al. 1996, 75. Schaefer fc Rin gwald 1995 . 76. Johnston Sz Kulkarni 1992 
78. 


! rampton, Cowley &; Fisher 1986, 79. Kraft 1964, 80. Mateo, Szkody &; Garnavich 1991 


69. 


995. 65. Shafter et al. 1995 


Casares et al. 1096 


horstensen 1997. 73. Howe ll. Li cbcrt fc Maso n 1994 


77. Rcinsch 1994, 


Spectral type uncertain since it was determined from a spectrum taken while the star was on the rise to outburst. 
Spectral type given as late G. 


(x) The width of the base of the emission lines. This gives _ lobe) and the mass of the primary star (e.g. Horne, Ver- 


the projected rotational velocity of the innermost parts of bunt & Schneider 1986). In combination with Kr, this can 


the accretion disc and hence a relation between Mi, Ri and 


i (e.g. Kuerster & Barwig 1988). This quantity is, however, 


difficult to measure with any precision. 

(xi) The separation of the two peaks in the emission lines, 
vd- This is used to relate the accretion disc radius (which 
must be less than the radius of the primary star’s Roche 


provide an upper limit to q. 

(xii) The orbital modulations of V sin i. With several as¬ 
sumptions about limb and gravity darkening these modu- 
lations can be mod e lled to constra in the inclination (e.g. 


Casares et al. 199f 


Sliahbaz 1998). This technique is re¬ 


ally a subset of mass derivations using Roche tomography 
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6 D. A. Smith and V. S. Dhillon 


Table 2. Spectral types, masses and radii of secondary stars in LMXBs. 


Star 

Porb ( h r ) 

Spectral Type 

M 1 (M q ) 

M 2 {Rq ) 

R 2 {Rq) 

References 

V518 Per (GRO J0422+32, Nova Per 1992) 

5.09 

M2+2V 




1,2 

MM Vel (Nova Vel 1993) 

6.86 

early K 




3 

V616 Mon (A0620-00) 

7.75 

K4V 

3.89-4.12 

0.19-0.32 

0.53-0.63 

4,5 

QZ Vul (GS 2000+25) 

8.26 

K3-6V 

6.04-13.9 

0.26-0.59 

0.62-0.81 

6 

GU Mus (Nova Mus 1991) 

10.4 

K3-4V 

6.98 + 1.45 

0.94 ± 0.40 

1.06 + 0.15 

7,8,9 

V2107 Oph (Nova Oph 1977) 

12.5 

K5V 




10 

V822 Cen (Cen X-4) 

15.1 

K5V 




11,12 

V1333 Aql (Aql X-l) 

19.0 

K5V 




11 

V1033 Sco (GRO J1655-40, Nova Sco 1994) 

62.9 

F3-6IV 

7.02 + 0.22 

2.34 + 0.12 

4.85 + 0.08 

13 

V404 Cyg (GS 2023+338) 

155.3 

K0IV 

12±| 

0 7 + 0-3 
u -' - 0.2 

6 0 +0 ' 7 
°’ U - 0.5 

14,15 


Clhevalier V Tlova.iskv 1 996 2. 

Casares et al. 1995. 3. Shahbaz et al. 19961. 4. McClintock &; Remillard 1986|. 

Marsh, Robinson & Wood 199 

l[ 6. Harlaftis. Horne &; FilinDenko 1996[ 7. |Orosz et al. 1996). 8. Casares et al. 1997 


9. This paper. 10. Harlaf tis e t al. 1997 
13. 


11 . 


Orosz & Bailyn 1997 14. Casares & Charles 1994, 15. 


Shahbaz, Naylor fe Charles 1996c[ 12. [vTcClintock fe Remillard 199C, 


Shahbaz et al. 1994b. 


(Rutten & Dhillon 1996), which utilise the modulations in 


mass determinations which are not listed in their catalogue. 
These are BT Mon, V895 Cen and V2051 Oph. We have also 
discovered a number of CVs in the literature which pos¬ 
sess measurements of Kr, Vsini, A (j>\/2 and i but which 
have had no mass calculations - for these objects (BD Pav, 
DX And, EX Dra, AM Her and the LMXB GU Mus), we per¬ 
formed Monte Carlo simulations similar to those described 


line strength and position as well as width to determine the 
component masses. 

(xiii) A mass-radius re l ation for the wh ite dwarf (e.g. 

Hamada & Salpeter 1961, Nauenberg 1972{ ). This is com¬ 
monly used in conjunction with item (vi). 

(xiv) A main-sequence mass-radius relation for the sec¬ 
ondary star. The assumption of a main sequence mass- 
radius relation is made frequently in the literature, and it is 
the validity of this assumption which we wish to test. 

Of these various measurements and techniques for mass 
determination, our preferred method is the combination of 
Kr, V sini and Acfri/2, or if the system is non-eclipsing, el¬ 
lipsoidal variations. These parameters are simple to measure 
if the secondary is sufficiently bright, do not depend on any 
assumptions, other than that the secondary fills its Roche 
lobe, and can have any biases (due to irradiation of the sec- 
onda ry, for example) corr ected for quite straightforwardly 
(e.g. |Davey fe Smith1991 ). 

Hitter & Kolb (1998) list values for the mass of the 
secondary star in 81 CVs. The vast majority of these mass 
determinations, however, have been derived using a main 
sequence mass-radius relation. As our goal is the derivation 
of a new mass-radius relation for CV secondary stars, we 
have been forced to omit from our consideration all of the 
mass estimates which make a main sequence assumption. 

This leaves only 20 mass determinations. Of these, 5 use 
the width of the emission lines (x and xi), which we believe 
to be extremely unreliable due to measurement uncertain¬ 
ties, 2 use the eclipse method (iv,v,vi) when it is clear from 
the shape of the light curve that this is invalid or at least 
unreliable, and 6 use determinations of Kw which are de- 
rived fr om radial velocity curves with significant phase shifts - secon dary stars ( Drew, Jones fc Woods 1993 , Casares et al 


by Smith, Dhillon & Marsh (1998) to calculate the system 
parameters - these are listed in Appendix C. 

In summary, we have 14 reliable CV mass determina¬ 
tions and 8 less reliable mass determinations, the latter de¬ 
pending on sometimes dubious Kw measurements, or else 
the uncertain assumption that the eclipse method is valid 
(UU Aqr). The masses and radii are listed in Table 3, and 
described in greater detail in Appendix A. Note also that a 
number of authors quoted only the masses of the secondary 
stars and not their radii - we have calculated the radii using 
equation (3). Where this is the case, the errors on the radii 
have been taken to be a third of the percentage errors on 
the mass (from equation 3). 

Th ere are just 5 LMXBs listed by Ritter & Kolb (1998, 
see also Beekman et al. 1997) which have reliable mass de¬ 
terminations. We have determined the system parameters of 
GU Mus using values of V sin i and Kr provided by Casares 
et al. (1997; who only determined upper and lower limits to 
the masses). The masses and radii of the LMXB secondary 
stars are listed in Table 2 and described in more detail in 
Appendix A. 

Those CVs with reliable mass determinations have been 
plotted as filled circles in Fig. 2, while those which use dubi¬ 
ous Kw determinations are ringed. The points representing 
DX And and AE Aqr are also ringed as they have evolved 


TO, 

eJu 


and fierce do not represent the motion of the white dwarf 
In fact, the only mass determination wholly dependent on 
a Kw me asurement whi ch we have not rejected is that of 
IP Peg, as Marsh (1988) was able to successfully correct for 
the phase shift in the radial velocity curve using the light 
centres method. 

The above filtering process leaves us w ith only 7 rcli- 
able CV mass determinations listed in the Ritter & Kolb 


199C). LMXBs are plotted as open triangles. Open circles 


rep resent the masses and radii of the 50 isolated stars listed 
by Andersen (1991) and Popper (1980), the former having 


(1998) catalogue. We have also uncovered a few additional 


smaller error bars than the latter. 

We performed least-squares power law and linear fits 
to the CVs using the formal error bars quoted in the lit¬ 
erature and listed in Table 3. We have only fitted to the 
unringed, filled circles in Fig. 2, i.e. DX And, AE Aqr and 
those systems which have dubious Kw measurements have 
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The secondary star in CVs and LMXBs 7 




Figure 2. The masses and radii of the secondary stars in CVs and LMXBs. The upper panel shows a power law fit (solid line) and a 
linear fit (dashed line) to the CV data (unringed, filled circles). The CV points which are ringed have been omitted from the fit. LMXBs 
are denoted by open triangles. Also plotted in the upper panel is the mass-radius relation (equation 7) derived by Warner (1995a,b). The 
lower panel shows the same data points as the upper panel, along with the theoretical models of Chabrier Sz Baraffe (1997, the dotted 
line) and the empirical relation derived by Clemens et, al. (1998. the solid li ne). The thick, solid line shows the secular evolution of the 
mass and radius of the secondary star computed by Kolb fc Baraffe (1998) . 50 stars in detached binaries with well-determined masses 
and radii have been plotted as open circles in the upper panel. See text for details. 


been omitted from the fits. The power law fit is plotted as 
a solid fine in the upper panel, and is given by 


R 

Rq 


(an ± cut,) (J^-) 


( 0 . 75 ± 0 . 04 ) 


(5) 


The linear fit is plotted as a dashed line in the upper panel, 
and is given by 

A = (0.93 ± 0.09) (^-) + (0.06 ± 0.03). (6) 


We have also plotted the mass-radius relation of Warner 
(1995a,b), 

R = M 13/15 , (7) 

in the upper panel of Fig. 2 as a dotted line, which is an 
approximate fit to the data set given by |Webbink (1990) . 
Of these three mass-radius relations (equations 5,6,7), we 
would recommend the use of our equation 5, which has been 
fitted to our data using a least-squares fitting procedure (as 
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8 D. A. Smith and V. S. Dhillon 


Table 3. Masses and radii of cataclysmic variables. Key: DN = dwarf nova; CN = classical nova; NL = nova-like; P = polar; IP = 
intermediate polar. Asterisks denote those stars whose system parameters were wholly dependent on Kyy measurements derived from 
radial velocity curves exhibiting significant phase shifts. 


Star Class P or b Mi q = M 2 /M 1 M 2 


R 2 Refs 


V2051 Oph 

DN 

1.50 

0.78 

± 

0.06 

0.19 

± 

0.03 

0.15 

± 

0.03 

0.16 

± 

0.01 

1 


OY Car 

DN 

1.51 

0.685 

± 

0.011 

0.102 

± 

0.003 

0.070 

± 

0.002 

0.127 

± 

0.002 

2 


EX Hya* 

IP 

1.64 

0.49 

± 

0.03 

0.19 

± 

0.03 

0.095 

± 

0.013 

0.15 

± 

0.01 

3, 4, 5 


HT Cas 

DN 

1.77 

0.61 

± 

0.04 

0.15 

± 

0.03 

0.09 

± 

0.02 

0.154 

± 

0.013 

6 


Z Cha 

DN 

1.79 

0.84 

± 

0.09 

0.20 

± 

0.02 

0.125 

± 

0.014 

0.172 

± 

0.010 

7, 8 


ST LMi 

P 

1.89 

0.76 

± 

0.30 

0.23 

± 

0.05 

0.17 

± 

0.07 

0.20 

± 

0.03 

5, 9, 10 

AM Her 

P 

3.09 

0.44 

± 

0.12 

0.64 

± 

0.10 

0.29 

± 

0.10 

0.33 

± 

0.04 

5, 11, 

12, 13 

IP Peg 

DN 

3.80 

1.15 

± 

0.10 

0.59 

± 

0.04 

0.67 

± 

0.08 

0.501 

± 

0.024 

14, 15 


CN Ori* 

DN 

3.92 

0.74 

r =t 

: 0.1 

0.66 

± 

0.04 

0.49 

± 

0.08 

0.46 

± 

0.03 

16 


UU Aqr* 

NL 

3.93 

0.67 

± 

0.14 

0.33 

± 

0.10 

0.20 

± 

0.07 

0.34 

± 

0.04 

17 


U Gem* 

DN 

4.25 

1.26 

± 

0.12 

0.46 

± 

0.03 

0.57 

± 

0.07 

0.510 

± 

0.023 

16, 18 


BD Pav* 

DN 

4.30 

0.95 

± 

0.10 

0.44 

± 

0.06 

0.43 

± 

0.10 

0.46 

± 

0.04 

5, 16, 

19 

DQ Her 

CN,IP 

4.65 

0.60 

± 

0.07 

0.66 

± 

0.04 

0.40 

± 

0.05 

0.49 

± 

0.02 

20, 21 


IX Vel 

NL 

4.65 

0.80+°;^ 

0.65 

± 

0.04 

0 52"*" 010 
u, ° —0.07 

0.530 

± 

0.025 

22 


V895 Cen 

P 

4.77 

0.93 

± 

0.17 

0.51 

± 

0.12 

0.48 

± 

0.17 

0.51 

± 

0.06 

23 


EX Dra 

DN 

5.04 

0.70 

± 

0.10 

0.84 

± 

0.12 

0.59 

± 

0.12 

0.59 

± 

0.04 

5, 24, 

25 

EM Cyg* 

DN 

6.98 

0.56 

± 

0.05 

1.35 

± 

0.16 

0.76 

± 

0.10 

0.79 

± 

0.04 

5, 26 


AC Cnc* 

NL 

7.21 

0.82 

± 

0.13 

1.24 

± 

0.08 

1.02 

± 

0.14 

0.92 

± 

0.05 

27, 28 


V363 Aur* 

NL 

7.71 

0.86 

± 

0.08 

0.89 

± 

0.03 

0.77 

± 

0.04 

0.83 

± 

0.04 

29 


BT Mon 

CN 

8.01 

1.04 

± 

0.06 

0.84 

± 

0.04 

0.87 

± 

0.06 

0.89 

± 

0.02 

30 


AE Aqr 

IP 

9.88 

0.79 

± 

0.16 

0.630 

± 

0.012 

0.50 

± 

0.10 

0.86 

± 

0.06 

5, 31, 

32 

DX And 

DN 

10.60 

0.51 

± 

0.12 

0.98 

± 

0.10 

0.50 

± 

0.14 

0.92 

± 

0.08 

5, 33, 

34 


References: 1. 

6 . 

10 . 

14. 

18. 

22 . 

25. 

28. 

30. 

33. 


Rantista, et, ah 1908 2. [Wood ct al. 198~. 3. Hollier 1096 , 4. ■it.erkeii et ak 1083 5. This na.ner 


CrooDcr 198f. 11. 

Davev & Smith 1996). 12. Southwell et al. 1995. 13. 

Wickramasinehe et al. 1991). 

Martin et al. 19££ 

15. Rarrera, V Vogt 19891 16. Friend et al. 1990a 

. 17. paDtista. Steiner &: Cieslinski 1994, 


Bcucrmann fc Thomas 199o|. 23 ~ Bucklcv ct al. 1998[ 24. |Billiiigton. Marsh .V Dhillon 1996 


Fiedler. RarwipjCJUan^ 26. gtover, Robin ^L^^atllgLT^^^-TL- ^lka.za.ki 


Schlcgcl. Kaitchuck fc Honeycutt 19sj. 29. Schlcgcl. Honeycutt fc Kaitchuck 1986 
Smith. Dhillon A: Marsh 1998| , 3 .. |Cas'arcs ct al. 1996 32. [Welsh, Horne &; Gomer 1993 


Drew, Jones &; Woods 199J, 34. Hilditch 1995 


Kitamura & Yamasaki 1982 


opposed to equation 7, where Warner fitted the Webbink’s 
data by eye, forcing the multiplicative constant to be unity 
and the power to be a simple ratio). 



stellar input physics for low-mass stars and brown dwarfs 
currently available. This particular sequence uses a white 
dwarf mass of Mi = 0.7Mq and an initial secondary star 
mass of M 2 = 0.6Mq. The main features in the upper part 
of the curve are the large radii at high masses (compared to 
the theoretical main sequence) due to thermal inequilibrium 
and the detachment of the secondary from its Roche lobe at 
M 2 = 0.22Mq, which is the point at which (in this model) 
the secondary becomes fully convective and magnetic brak¬ 
ing is assumed to cease. During this detached phase, there 
is a 25 per cent fall in radius as the secondary star relaxes 
back to thermal equilibrium. 


The lower part of the curve follows the theoretical main 
sequence, with the secondary star close to thermal equilib¬ 
rium at low mass transfer rates, until the minimum mass for 
nuclear fusion is reached and the secondary becomes a de¬ 
generate, brown dwarf-like object with R 2 oc Mj -1 ^ 3 . Note 
that the evolutionary sequence should be used with caution 
when interpreting the observed CV masses, as it represents 
the change in mass and radius of a single secondary star 
with time, whereas the observed data is a snapshot of the 
masses and radii of the population of secondary stars at this 
moment in time. To make a proper comparison with theory, 
one would need to compute the evolutionary tracks of a sam¬ 
ple of secondary stars generated by a population synthesis 
code. 


As a group, the CVs appear to lie slightly above the the¬ 
oretical main sequence of Chabrier & Baraffe (1997) but ap¬ 
pear to fit the main sequence as defined by the distribution 
of detached systems very well. There are a few departures; 
as already noted, DX And and AE Aqr have evolved sec¬ 
ondary stars and lie above the main sequence, while UU Aqr 
and IP Peg both lie above the main sequence, being under¬ 
sized for their masses, but within the error bars they are 
still consistent with the main sequence. Rather surprisingly, 
those systems which we have omitted from the fitting pro- 
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The secondary star in CVs and LMXBs 9 


cess because of dubious Kw measurements also follow the 
main sequence very closely. Unfortunately, the errors on the 
masses and the scatter of the points in Fig. 2 do not allow 
us to say if the secondary stars in CVs follow the disrupted 
magnetic braking model, or the kinked lower main sequence 


Clemens et al. (1998)| .Plrhis means our data are unable 


of 

to test which of these two period gap formation mechanisms 
is correct (see Section 7). There is also no evidence in these 
data for systems which have evolved beyond the orbital pe¬ 
riod minimum and h ave low-mass, degenerate brown dwa rf- 
like secondaries (see Howell, Rappaport & Politano 1997). 

Of the three short period LMXBs, V616 Mon and 
QZ Vul both contain secondary stars which lie above the 
main sequence and are somewhat evolved, while GU Mus 
lies just on the main sequence (to within the error). The 
secondary stars in the long period LMXBs V404 Cyg and 
V1033 Sco are both evolved (V1033 Sco is also very massive) 
and do not appear in Fig. 2. 


4 THE MASS-SPECTRAL TYPE RELATION 


An assumption often made in the estimation of system pa¬ 
rameters in CVs is that the spectral type of the secondary 
star can be used to estimate its mass. In Fig. 3 we have plot¬ 
ted those CVs and LMXBs which have measured masses and 
spectral types given in Tables 1, 2 and 3. Also plotted are the 


detached systems from Andersen (1991) and Popper (1980) 


The mass of the secondary is thought to be the most im¬ 
portant factor in determ ining its effective temperature, and 
therefore spect ral type (King & Kolb 1998). Stelile, Ritter 
& Kolb (1996) show that simple homology relations lead to 
a relationship between the effective temperature of the star 
(and therefore spectral type) and the radius of the form 
T e ff oc P'', where |ju| « 1. Stars should therefore expand 
or contract without changing their spectral types by much 
and the spectral type of a star should be almost completely 
dependent on its mass. 

Fig. 3 shows that this is not the case, and that there 
is a huge range in mass for a given spectral type. Around 
M5, for example, lie IP Peg (M4.5) and Z Cha (M5.5); the 
secondary star in IP Peg is five times more massive than 
that in Z Cha. Among the LMXBs there is also little or 
no correlation between spectral type and masses; GU Mus 
(K3-4V) is around four times more massive than V616 Mon 
(K4V). The large scatter in masses for a given spectral type 
are, however, shared by the detached systems (e.g. around 
M4-5 and K0). So although the secondary stars in CVs and 
LMXBs are not too different from main sequence stars, there 
is too much variation to use the spectral types of CV and 
LMXB secondaries to estimate the secondary star masses, 
as is also the case with the detached stars. 


T This has recently been shown not to create a period gap at all, 
but rather two spikes at P ~ 2 hr and P ~ 3 hr with the prob¬ 
ability of discovering a CV in the van no lower than discovering 
one outside the gap (Kolb, King fe Ritter 1998|). 



Spectral Type 

Figure 3. The masses and spectral types of the secondary stars 
in CVs and LMXBs. The filled circles represent the CVs, the tri¬ 
angles the LMXBs, and the unfilled circles the detached systems. 


5 THE MASS-PERIOD AND 

RADIUS-PERIOD RELATIONS 

Another relation of interest is that between the orbital pe¬ 
riod and the mass of the secondary star. We have used the 
masses and periods listed in Table 3 to derive a mass-period 
relation. By fitting a power law to the data (excluding the 
evolved secondaries in DX And and AE Aqr and the sys¬ 
tems with dubious Kw determinations), we have derived 
the following relationship: 

= (0.038 ± 0.003)P (1 ' 58±o '° 9) . (8) 

We have also performed a linear fit and derived the following 
mass-period relationship, 

M 

— = (0.126 ± 0.011)P - (0.11 ± 0.04). (9) 

Warner (1995a,b) derived the following semi-empirical 
mass-period relation, which utilised his mass-radius rela¬ 
tion (equation 7): 

-^- = 0.065P 5/4 . (10) 

Mq 

The mass-period relation is shown in the left-hand 
panel of Fig. 4. The filled circles represent the CVs; those 
which are ringed have been omitted from the ht either be¬ 
cause they depend upon dubious Kw measurements, or be¬ 
cause they are evolved to some degree. The solid line is the 
power law ht (equation 8), the dashed line is the linear ht 
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Orbital Period (hr) Orbital Period (hr) 

Figure 4. The masses and radii of the secondary stars in CVs versus their orbital periods. In each panel, the solid line is the power law 
fit to the CV data (equations 8 and 11) and the dashed line is a linear fit to the CV data (equations 9 and 12). The dot ted lines are the 
semi- emp irical relations of Warner (1995a,b; equations 10 and 13). The thick solid line is the evolutionary sequence of Kolb & Baraffe 
(1998)|. [The ringed points have been omitted from the fits - see text for details. 


(equation 9), and the dotted line is W arner’s relation (equa¬ 


tion 10). The evolutionary sequence of Kolb &; Baraffe (1998) 
is plotted as a thick, solid line. 


-P- = (0.081 ± 0.019)P (118±0 04) . (11) 

it© 

We have also performed a linear fit, which is given by 


The linear fit (unweighted rms deviation = 0.10) is sur¬ 
prisingly better than the power law fit (rms = 0.12); one 
would have expected a homologous relation between mass 
and period, and therefore a power law to be the best fit. 
Warner’s fit is also surprisingly good (rms = 0.11), but it 
is a poorer fit to the systems with short periods and those 
with the most accurate mass determinations. The evolved 
systems DX And and AE Aqr lie a long way from the fit, 
while IP Peg is twice as massive as the fits predict. The 
LMXBs all have lower masses than the CV fits predict. The 
period gap is represented in Fig. 4 by the flat section of the 
evolutionary sequence at M = 0.22M@; the data are insuffi¬ 
ciently accurate or numerous to say whether the secondary 
stars follow it. There are no stars lying on the degenerate sec¬ 
ondary star arm of the evolutionary sequence implying that 
none of the short period systems in our sample containing 
the degenerate, brown dwarf-like secondaries we expect in 
post-period minimum CVs. 

A similar procedure to that above has been followed to 
derive a radius-period relation for CV secondary stars. A 
power law fit has been applied to the data in Table 3, again 
excluding DX And and AE Aqr and the systems with dubi¬ 
ous Kw measurements. The resulting radius-period relation 
is 


-p- = (0.117 ± 0.004)P - (0.041 ± 0.018). (12) 

it© 

The corresponding semi-empirical radius-period relation of 
Warner(1995a,b) is 


-P- = 0.094P 13/12 . (13) 

ft© 

The right hand panel in Fig. 4 shows the radius-period 
relation for CVs. The solid line is the power law fit (equa¬ 
tion 11) to the CV data (unringed circles), the dashed line is 
the linear fit (equation 12) and the dotted line is Warner’s re¬ 
lation (equa tion 13). The thick solid line is the evolutionary 
sequence of Kolb & Baraffe (1998), The rms scatters from 
the fits show them to be of similar quality; 0.033 for the lin¬ 
ear fit, 0.042 for the power law and 0.036 for Warner’s rela¬ 
tion. The correlation here is stronger than that in the mass- 
period diagram, as is to be expected: from equation (3), any 
deviation in mass from the mass-period relation must be 
matched by a deviation one third the size in the radius- 
period diagram. The LMXBs have smaller radii than their 
orbital periods predict, again the deviation from the fit is 
around one third of that in the mass-period plot. 
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Figure 5. The masses of the white dwarfs and the mass ratios 
of CVs plotted against orbital period. Open symbols are used for 
those mass determinations which used unreliable Kyy measure¬ 
ments. 


6 WHITE DWARF MASSES 

The careful selection of systems with well measured sec¬ 
ondary star masses and radii also provides us with the most 
accurately determined sample of CV white dwarf masses. In 
Table 4, we list the mean white dwarf mass for each of the 
CV subtypes, above and below the period gap, along with 
the standard errors, followed by the number of systems in 
parentheses. We have computed weighted and unweighted 
means for the white dwarfs - the methods used are given 
in the footnotes to Table 4. The unweighted mean is prob¬ 
ably the more reliable statistic to use, since the range of 
formal errors in the white dwarf mass determinations is so 
great; the photometric determinations which have small for¬ 
mal errors but certainly much larger systematic errors (e.g. 
OY Car) dominate in the calculation of the weighted means. 
The dispersions on the unweighted means are also more rep¬ 
resentative than those quoted for the weighted means. White 
dwarf masses are plotted against orbital period in the up¬ 
per panel of Fig. 5, with each CV subtype represented by a 
different symbol. 

The systems below the period gap have a lower mean 
mass than those above, which is to be expected, since high 
mass white dwarfs are required to support stable conserva¬ 
tive mass transfer from the more massive secondaries which 
reside in long period CVs. Most of the white dwarf masses 
are consistent with CO white dwarfs; the white dwarfs in 
U Gem and IP Peg are near the minimum mass for ONeMg 


white dwarfs (~ 1.15Mg, e.g. Iben, Ritossa & Garcia-Berro 
1997), but none are light enough to be He white dwarfs 


(<0.4Mq). The mean mass for systems below the gap is 
probably more representative of the class as a whole, since 
strong selection effects come into play: it is thought that up 
t o 99 per ce nt of the actua l CV population are below the gap 
(Kolb 1993 - but see also Patterson 1998), despite the fact 
that they comprise only about 40 per cent of those observed. 
The white dwarf mas ses show no trends with CV sub-type. 

Web bulk (1990) performed statistics on a larger sample 
of white dwarf masses, but a large number of these deter¬ 
minations were unreliable, depending on measurements of 
emission line profiles (FWHM, double-peak separation, rms 
line-widths) which had been calibrated against a few well 
observed double-lined CVs and Algols. It is, however, worth 
comparing Webbink’s figures with ours as his mean white 
dwarf mass es are the most o ften quoted. The mean masses 
derived by Webbink (1990) for all systems, and for those 
above and below the period gap, are listed in Table 4. We 
reach the same general conclusions as Webbink, the mean 
white dwarf mass above the period gap is higher than that 
below, and that the overall mean mass is 0.76 ± 0.22A/q, 
hig her than the mean mass of single wh ite dwarfs (0.6 Mq, 
e.g. Bergeron, Liebert fc Fulbright 1995| ). 

The mass ratio-orbital period diagram in the lower 
panel of Fig. 5 shows the near linear relation we expect, 
since the secondary stars’ masses are well correlated with 
their orbital periods. The scatter in this is largely due to 
the scatter in white dwarf masses. Again there is no trend 
according to subtype; the two systems with anomalously 
large mass ratios, EM Cyg and AC Cnc, are a dwarf nova 
and a novalike respectively. Note that this sample of accu¬ 
rate mass ratio determinations is not exhaustive, there exist 
several systems which have the mass ratio determined, but 
which do not have accurate mass determinations (e.g. be¬ 
cause the inclination is unknown.) 

In our sample, the unweighted mean value of the mass 
ratio (performed on those systems which do not have du¬ 
bious Kw values) is q = 0.17 ± 0.05 below the period gap 
and q = 0.70 ± 0.15 above the period gap. These values 
are consi stent with the statistical analysis of eclipse dura¬ 
tions by Bailey (1990), who found that the mean value of 
the mass ratio should be q = 0.13 ± 0.03 below the gap and 
q = 0.65 ± 0.12 above the gap. 


7 DISCUSSION 

7.1 Are CV secondaries main sequence stars? 

The spectral type-period diagram (Fig. 1) shows that the 
secondary stars in CVs lie very close to the main sequence 
defined by the detached systems. There is little differerence 
between the distribution of CV secondaries and that of the 
detached stars, up to an orbital period of about 7-8 hours. 
Beyond that, for the reasons stated in section 1, the sec¬ 
ondary stars are required to be evolved, and we see AE Aqr 
(P = 9.88 hr) and DX And (P = 10.60 hr) lying well above 
the rest of the distribution with later spectral type than the 
fit predicts. 

The mass-radius diagram (Fig. 2) clearly shows that 
the secondary stars in CVs are on the whole indistinguish- 
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Table 4. The mean white dwarf masses for each of the CV subtypes, for systems below the period gap, above the period gap, and in 
total. The standard errors are also listed, with the numbers of systems used in determining the mean given in parentheses. Note that 
DQ Her has been counted as both a classical nova and an intermediate polar. 


A. Unweighted Average 0, 

Sample 

Below period gap 

Above period gap 

All periods 

Dwarf novae 

Nova-likes 

Polars 

Intermediate polars 
Classical novae 

All systems 

0.73 ±0.10 (4) 

0.76* (1) 

0.49* (1) 

0.69 ±0.13 (6) 

0.84 ±0.29 (7) 

0.79 ± 0.08 (4) 

0.69 ± 0.35 (2) 

0.70 ±0.13 (2) 

0.82 ±0.31 (2) 

0.80 ±0.22 (16) 

0.80 ±0.24 (11) 
0.79 ±0.08 (4) 
0.71 ±0.25 (3) 
0.63 ±0.15 (3) 
0.82 ±0.31 (2) 
0.77 ±0.21 (22) 

B. Weighted Average b 

Sample 

Below period gap 

Above period gap 

All periods 

Dwarf novae 

Nova-likes 

Polars 

Intermediate polars 
Classical novae 

All systems 

0.68 ±0.01 (4) 

0.76 ±0.30 (1) 
0.49 ±0.03 (1) 

0.66 ±0.01 (6) 

0.75 ±0.03 (7) 

0.82 ± 0.06 (4) 

0.58 ± 0.09 (2) 

0.63 ± 0.06 (2) 

0.85 ±0.05 (2) 

0.78 ±0.02 (16) 

0.69 ±0.01 (11) 
0.82 ± 0.06 (4) 
0.60 ±0.09 (3) 
0.52 ± 0.03 (3) 
0.85 ± 0.05 (2) 
0.68 ±0.01 (22) 

Webbink (1990) 

All systems 

0.66 ±0.01 (26) 

0.81 ± 0.04 (58) 

0.74 ±0.04 (84) 


*No associated error. 


a M 
b M 



E N Mj /sr^N i 
i =1 cr? ' 2-^i= 1 a-? 


J M 


t 2 
’ M 




able from the observed main sequence stars in detached bi¬ 
nary systems. Again, the only outlying points are AE Aqr 
and DX And, which almost certainly have evolved sec¬ 
ondary stars. The most interesting question from a mass- 
determination point of view, however, is not whether sec¬ 
ondary stars as a whole are main sequence stars, but rather 
at what point can we no longer apply a main sequence mass- 
radius relationship. From Figs. 1 & 2 we can say that CVs 
with periods of up to 7-8 hours almost always have main 
sequence secondary stars, although one must always beware 
of the existence of systems with peculiar secondaries and 
use the mass-radius relation with caution (e.g. by correctly 
propagating the errors in our mass-radius relation when de¬ 
termining masses). 


We have seen that the expected divergence from the 
main sequence does occur in long periods CVs, but there 
is as yet no firm evidence for post-period minimum CVs 
with degenerate brown dwarf-like secondaries in short period 


CVs. The evolutionary sequence of Kolb & Baraffe (1998) 


shows us where to look, but as yet there have been no brown 
dwarf-like secondaries detected in CVs either by their spec¬ 
tral types (e.g. by looking for spectral types much later than 
predicted by equation 4) or by their position on the mass - 
radius diagram (but see Howell, Hauschildt & Dhillon 1998). 
This does not mean, of course, that there are no post-period 
minimum CVs; these systems should be intrinsically very 
faint, and near the period minimum (where WZ Sge and 
many other suspected post-period minimum CVs lie) their 
secondary stars will be difficult to distinguish from late-type 
main sequence stars. 


7.2 Are LMXB secondaries main sequence stars? 


The spectral types of the LMXBs seem to follow the same 
distribution as the CVs with orbital period, while their radii 
are larger than main sequence stars of the same type, sug- 
gesting they are eit her somewhat evolved (as predicted by 
King & Kolb 1997), or else have expanded towards a new 


state of thermal equilibrium due t o the level of X-ray i rradi- 
ation to which they are exposed (Podsiadlowski 1991). The 
upper orbital period limit above which secondary stars can¬ 
not be on the empirical main sequence rises (from around 
9 hours for CVs), to around 15 hours, using the theoreti¬ 
cal upper neutron star mass limit of ~ 3Mg with a typical 
AO main seq uence star of mass 2.40Mg and radius 1.87Mg 
( |Gray 1992 ). In LMXB black hole candidates, of course, 


there is no such limit to the mass of the black hole and 
therefore no limit on the secondary star. 


7.3 Consequences for the disrupted magnetic 
braking model 

The generally accepted mechanism for the formation of the 
period gap in CVs is the disrupted magnetic braking model. 
Above the period gap magnetic braking is the dominant 
mechanism for the angular momentum loss which drives 
mass transfer. As the CV moves towards a period of 3 hours, 
the secondary star’s mass falls to around 0.25Af,7i. at which 
point it becomes (almost) fully convective ( Kolb 1993 ). Mag¬ 
netic braking is then thought to cease, allowing the sec¬ 
ondary, which due to its high mass transfer rate has been 
out of thermal equilibrium, to relax and shrink inside its 
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Roche lobe. This cuts off mass transfer (and the cataclysmic 
behaviour) until further loss of angular momentum via grav¬ 
itational radiation reduces the period to around 2 hours, at 
which point the secondary comes back into contact with its 
Roc he lobe and recommences mass tran sfer (Spruit & Ritter 
1983; Rappaport, Joss & Verbunt 1983). 

A mass-radius diagram would be expected to show a 
“kink” at around a mass of 0.25Mq, as the CV secondary 
stars revert from being out of thermal equilibrium and there¬ 
fore undermassive for their radii above the period gap, to be¬ 
ing thermally relaxed at low masses below the period gap. 
This is shown in the secular evolution model plotted in the 
lower panel of Fig. 2. If a star enters the period gap at 
Pabove = 3.0 hr and emerges at Pbelow = 2.0 hr, without 
having lost any further mass, then, using equation (3), its 
radius must have shrunk by a factor 


( Pbelow \ _ / Pbelow \ ^ _ q 

Pabove / V Pabove / 


(14) 


i.e. a decrease of almost 25 per cent. A narrower period gap, 
e.g. 2.2-2.8 hrs, would see a smaller drop in the radius at 
the top of the period gap, in this case 15 per cent. 

We are unable to see the expected sudden drop in Fig. 2, 
although the data is sparse in the crucial region around 
0.25 Mg, the only system there being AM Her. In this in¬ 
teresting region, just above the period gap, there is a dearth 
of dwarf novae, instead the CV population is dominated b y 
peculiar nova-likes such as the SW Sex stars ( Dhillon 1996 ), 
which have intrinsically bright discs and almost undetectable 
secondary stars, and for which determination of system pa¬ 
rameters is therefore very difficult. 

Clemens et al. (1998) see a dip in the colour-magnitude 
diagram for M-dwarfs at My ~ 12 (around 0.25Mq). They 
suggest that this means that the radius of the secondary 
star in a CV would shrink more rapidly in this region for 
the same rate of mass loss, leading to the rapid crossing 
of the 2-3 hr orbital period range. This would then be a 
possible mechanism for the formation of the period gap. This 
interpr etation has been disputed by Kolb, King & Ritter 


(1998), who show that if the mass-radius relation of Clemens 
et al. held, then rather than producing a period gap, the 
period distribution would have two spikes at the upper and 
lower edges of the “gap”, with the probability of discovering 
CVs inside the “gap” the same as discovering them outside 
the “gap.” 

Unfortunately we are not yet at the stage where we can 
confirm or disprove the disrupted magnetic braking model 
observationally. Our dataset is simply too sparse and insuffi¬ 
ciently accurate. Several more accurate measurements of sec¬ 
ondary star masses and radii are needed, especially around 
the period gap (in which there are sadly few secondary star 
detections) in order to differentiate between these competing 
period-gap formation theories. 


8 CONCLUSIONS 

(i) We find there are a total of 55 reliable spectral type 
determinations and only 14 reliable mass determinations of 
CV secondary stars (10 and 5, respectively, in the case of 
LMXBs). 

(ii) We have derived new mass-radius, mass-period, 


radius-period and spectral type-period relations for CV sec¬ 
ondary stars, using a carefully selected sample of CVs with 
well measured system parameters. 

(iii) The secondary stars in CVs with periods below 7-8 
hours are, as a group, indistinguishable from main sequence 
stars in detached systems in terms of spectral type, mass 
and radius. 

(iv) The secondary stars in LMXBs show some evidence 
for evolution, with radii which are slightly too large for their 
masses. 

(v) We have shown that the assumption that the spectral 
type of the secondary star in CVs and LMXBs provides a 
good estimate of its mass is not a good one. 

(vi) We have calculated the mean white dwarf mass in 
CVs, for the various CV subtypes, both above the period 
gap (where we find M\ = 0.80 ± 0.22M 0 ) and below the 
period gap (where we find M\ = 0.69 ± O.I 3 M 0 ). 

(vii) We have shown that accurate measurements of 
masses and radii in CV secondary stars can be used to con¬ 
strain CV evolution and provide evidence for or against the 
disrupted magnetic braking theory. 

(viii) We have demonstrated the need for many more ac¬ 
curate CV mass measurements, especially around the pe¬ 
riod gap, to test the disrupted magnetic braking model, and 
among the faint short period CVs to search for the predicted 
post-period minimum systems. 
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APPENDIX A 


Notes on individual systems 


(i) V2051 Oph. Baptista et al. (1998) used a purely pho¬ 
tometric method to find the masses of the components in 
this dwarf nova. The white dwarf and bright spot eclipses 
are both clearly visible in the light curve so the additional 
measurement of Ac/q /2 allowed the system parameters to be 
determined using the whit e dwarf mass-radius relation of 
Hamada fc Salpeter ( 1961) and a Monte Carlo simulation. 

(ii) OY Car. Wood et al. (1989) also used the photomet¬ 
ric method to determine the masses in this dwarf nova. Sev¬ 


eral other mass determinations cited by Wood et al. (1989) 
all give higher values for M 2 . 

(iii) EX Hya. This interm ediate polar has a detectable 


secondary star. Hellier (1996) combined Kw (Hellier et al 


1987) and Kr (measured using skew mapping; smith 


Cameron & Tucknott 1993) with A(f)ii 2 to obtain the mass. 


We have calculated the radius of the secondary star using 
equation (3). 


(iv) HT Cas. Horne, Wood fe Steining (1991) performed 
a similar analysis to that of Wood_et al. (1989)|. 


(v) Z Cha. Wood et al. (1986) 


used the same tech¬ 
nique as that applied to V2051 Oph, and obtained values 
of 0.081 ± 0.003Mq and 0.149 ± 0.0027?q for the mass and 
radius of the secondary star respectively. Wade & Horne 


(1988) derived another mass estimate, using the mass ratio 


and inclination i derived by Wood et al. (1986) , but rather 
than using the white dwarf ingress/egress timings, they ob- 
tained a measurement of Kr. The higher values derived by 
Wade fe Horne (1988) are those used throughout this pa- 
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per. The inconsistency between the two techniques cannot 
be satisfactorily explained by errors in the Hamada-Salpeter 
relation or the presence of a thick boundary layer above the 
surface of the white dwarf. 


(vi) ST LMi. Shahbaz & Wood (1996) measured the sec¬ 
ondary star’s radial velocity in this non-eclipsing polar, and 
determined the masses us ing the value o f i derived from po- 
larimetric measurements (Cropper 1986). 


(vii) AM Her. We have calculated the mass and radius 
of the secondary in the prototypic al polar, using the values 
of I\ n and V sin i determined by 5outhweU_gt_al_(1995)^ 
with the appropriate K-correction given by pavey fe Smith 
£1996^, and i determined using polarimetry by Wickramas- 
inghe et al. (199~j~j| . Southwell et al. (1995) found an 8 per 
cent bias in their measurement of V sin *, we therefore adopt 
V sin* = 100±10 kms -1 . Note that Shahbaz & Wood (1996) 

-l 


found V sin i = 68 ± 12 km s~ , using fewer spectra, and did 
not correct for irradiation effects. We therefore omit this 
measurement from the calculation. See Appendix C. 

(viii) IP Peg. 


Maxsh£1988^ applied the light centres 
method to obtain Kw- [Martin et al. ( 1989)| found Kr and 
A(f> 1/2 which provides a direct spectroscopic determination 
of the mass and r adius of the secondary star. 

(ix) CN Ori. Friend et al. (1990a)| measured Kr, and 
combined this with possibly unreliable values of * and Kw 


derived from a radial velocity study by Mantel et al. (1987) 
the details of which are unpublished. _ 


(x) UU Aqr. Baptista, Steiner fe Cieslinski (1994) used 
the method of Wood et al. (1989) to calculate the system 


parameters of this nova-like from photometry of the eclipses. 
Unfortunately the eclipse light curves do not have the shape 
usually required for the eclipse technique to be completely 
valid. We therefore group this mass determination with the 
possibly unreliable Kw determinations. 

(xi) U Gem. The value of Kr found by Friend et al 


(1990a) was combine d with Kw m easured by Stover (1981a) 


and i estimated by Smak (1976) 


from the location of the 


bright spot to calculate the masses. 


(xii) BD Pav. Kr and V sin * were measured by Friend 


et al. (1990a) 


and have been combined with the measure¬ 


ment of Acpi /2 ( Harrop-Allin fe Warner 1996 and references 
therein) to obtain the mass and radius of the secondary star. 
See Appendix C. 

(xiii) DQ Her. Horne, Welsh & Wade (1993) measured 
Kr and V sin* and combined these with previous measure¬ 
ments of Kw and A <f>i /2 to obtain a full set of system pa¬ 
rameters, using a Monte Carlo simulation. 

(xiv) IX Vel. Beuermann & Thomas (1990) detected 


emission lines emanating from the secondary star in this 
bright, non-eclipsing nova-like, which effectively make it a 
double-lined binary. Using kinematic and geometric mod¬ 
elling, and using the Balmer line light curves to constrain *, 
they obtained the mass and radius of the secondary star. 

(xv) Vi895 Cen (=EUVE J1429 + 38.0). The mass and 
radius of the seconda ry have recently been determined by 
Buckley et al. (1998) . Kr and V sin * were obtained using 
the Nal 8190A absorption doublet and i was found from 
ellipsoidal variations. 


(xvi) EX D r a (= HS 1804 + 6753). Fiedler, Barwig & 
Mantel (1997) measured Kw and Kr and obtained * from 
the geometry of the white dwarf and bright spot eclipses. 
However the radial velocity curve of the Ha emission line is 


out of phase by ~ 0.2 and the value of Kw determined is 
therefore unreliable. We have recomputed the system param¬ 
eters using a Monte Ca rlo simulation with the values of Kr 
and V sin * derived by jBillington, Marsh fc Dhillon (1996) 


and A^i /2 derived by Fiedler, Barwig fe Mantel (1997)| . See 
Appendix C. 

(xvii) EM Cyg. This double-lined dwarf nova has glanc¬ 
ing eclipses, and ha s had the radial velocities of both co mpo- 
nents measured by Stover, Robinson & Nather (1981). Un¬ 
usually, the secondary star is more massive than the white 
dwarf. This means that EM Cyg lies outside the ( M 2 ,q ) 
range of thermal and dynamical stability of mass transfer, 
assuming standard properties of Population I stars. These 
are, however, sensitive to metallicity, opacities and convec¬ 
tion theory used to compute the models. The fact that 
EM Cyg is observed to be in a moderate state of mass trans¬ 
fer means the star is probably not unstable. 

(xviii) AC Cnc. Another eclipsing double-lined system 
(^chlegel, Kaitchuck fc Honeycutt 1981), also with q > 1. 


(xix) V363 Au r (Lanning 10). Schlegel, Honeycutt & 


Kaitchuck (1986) obtained the mass and radius of the sec¬ 
ondary star through measurements of Kr and Kw and es¬ 
timated i using A(f> 1 / 2 . _ 


(xx) BT Mon. Smith, Dhillon & Marsh (1998) 


measured 

Kr from the weak secondary star absorption lines using 
skew mapping, Kw, V sin * and A(f > 1 /2 to obtain a full set 
of system parameters using a Monte Carlo simulation. 


(xxi) AE Aqr. Casares et ah (1996) detected absorption 
features from the secondary star and obtained measurements 
of Kr and V sin i. By modelling the way V sin * changes with 
orbital phase, they were able to constrain i. 

measured 
used 


(xxii) DX And. Drew, Jones fe Woo ds_(1993) 

Kr and V sin i and also estimated Kw - |Hilditch (199~5) 


ellipsoidal variations to provide an estimate of i. We have 
used this estimate with the measurements of Drew et al. 
to calculate the mass and radius of the secondary star. We 
pre fer the q determinat ion of Drew et al. over that derived 
by Bruch et al. (1997) because Bruch et al. used a phase- 
shifted Kw measurement from poorly wavelength calibrated 
data. 


(xxii i) V616 Mon (A0620-00). Marsh, Robinson & Wood 


(1994) measured Kr, V sin* and ellipsoidal variations in the 


equivalent width of the Ha emission line. We have used 
their la errors. A further constraint on * is given by the 


grazing eclipses observed by Haswell et al. (1993) , giving 
i ~ 70°. However, i nfra-red photometry by Shahbaz, Naylor 
fc Charles (1994a) of the ellipsoidal variations give a lower 
inclination, i ~ 40°, which implies a higher mass and radius, 
M 2 ~ 0.6M© and R 2 ~ 0.8Rq. 

(xxiv) QZ Vul (GS 2000+25). Harlaftis, Horne & Filip- 


penko (1996)| measured V sin i and Kr and combined these 
measurements with a wide range of allowable values of * to 
obtain an estimate of the mass of the secondary. 

(xxv) GU Mus (Nova Mus 1991, GRS 1124-68). Casares 


et al. (1997) measured Kr and V sin* from the secondary 
absorption lines. * has been e stimated from the ellipsoidal 
variations ( Orosz et al. 199^ ). Using these parameters we 
have calculated the mass and radius of the secondary using 
a Monte Carlo simulation. See Appendix C. 

(xxvi) V1033 Sco ( Nova Sco 1994, GRO J1655-040). 
Orosz fe Bailyn (1997) measured Kr spectroscopically and 
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modelled the light curves, which contained substantial ellip¬ 
soidal variations, to obtain q and i. _ 

(xxv ii) V 404 Cyg (GS 2023+338) Casares fc Charles 


£1994^ measured V sin i and Kr spectroscopically. Shahbaz 


et al. (1994b) modelled the ellipsoidal variations to obtain 


i and hence the system parameters. See also Shahbaz et al 


(1996a) 


APPENDIX B 


Simple formulae for estimating the mass ratio and 
secondary star mass 


While the best method of calculating the masses of primary 
and secondary stars and the other system par ameters in CVs 
and LMXBs is the Monte Carlo method (e.g. Smith, Dhillon 
& Marsh 1998), simple approximations can give fairly pre¬ 
cise estimates of q and M 2 . 

Kr and V sin i are respectively given by 


Kr = — -———- sin i 

P (1 + 9) 

(15) 

and 


T/ • • 271 - . . 

V sin 1 — —TC2 sin 1. 

(16) 

Combining these two equations gives 


R 2 ,, , Fsini 

a (1 + s)= K R ' 

(17) 


which using equation (2) yields a simple cubic formula for 
approximating q, 


q( 1 + q) 2 = 9.6 


( 


Vsini\ 

Kr ) 


3 


(18) 


Also, the density-period relation given by equation (3) 
can be used in combination with equation (16) to give a 
simple estimate of the mass of the secondary in terms of 
V sin* and P if i is known roughly 



0.042 ( 


Psini \ 3 P(hr) 
100 km s -1 / sin 3 i 


(19) 


For eclipsing systems, sin 3 i can be approximated as 
0.98 without introducing a significant error. This then leads 
to the mass of the secondary being solely dependent on 
P sin i and P. 



0.043 ( 


P sin i 
100 kms -1 



( 20 ) 


APPENDIX C 

Results from Monte Carlo simulations 

A number of the secondary star masses and radii have been 
derived using Monte Carlo simulations with values of P, 
Kr, V sini, i and A</) 1 / 2 recovered from the 
similar to that described by 
. The resulting system parameters are given 


method used is 
fe Marsh (1998) 
in Table 5. 


literature. The 
Smith, Dhillon 
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Table 5. Results from the Monte Carlo simulations. Parameters in bold face were those input to the Monte Carlo simulator. 



AM Her 

BD Pav 

EX Dra 

DX And 

GU Mus 

P (hr) 

3.09 

4.30 

5.04 

10.60 

10.38 

K r (kms -1 ) 

179±2 

278±4 

210±14 

105.8±3.8 

420.8±6.3 

K w (kms -1 ) 

115 ± 18 

123 ± 17 

176 ± 19 

103 ±9 

56 ± 17 

V sin* (kms -1 ) 

100±10 

125±10 

140±10 

79±5 

106±13 

i 

52±5 

73.4 ± 0.9 

82.1 ± 2.0 

49±4 

59.5±5.5 

A</>i/ 2 

- 

0.040±0.006 

0.1103±0.0001 

- 

- 

Mi (M©) 

0.44 ± 0.11 

0.95 ±0.10 

0.70 ±0.10 

0.51 ±0.12 

6.98 ± 1.45 

m 2 (m 0 ) 

0.29 ±0.10 

0.43 ±0.10 

0.59 ±0.12 

0.50 ±0.14 

0.94 ± 0.40 

R 2 (-R©) 

0.33 ±0.04 

0.46 ± 0.04 

0.59 ± 0.04 

0.92 ±0.08 

1.06 ±0.15 

<? 

0.64 ±0.10 

0.44 ± 0.06 

0.85 ± 0.12 

0.98 ±0.10 

0.13 ±0.04 

a ( Rq ) 

0.96 ±0.09 

1.48 ±0.05 

1.58 ±0.08 

2.44 ±0.19 

4.77 ±0.34 


© 0000 RAS, MNRAS 000 , 000-000 






